Geometrical locus of massive test particle orbits in the space 
of physical parameters in Kerr space-time. 

O ! F. Fayos*, Ch. Teijon 

^ ! Department of Applied Physics, UPC, Barcelona, Spain. 

I— J ; 

I ■ 

Abstract 

Gravitational radiation of binary systems can be studied by using the adiabatic approxima- 
^yi' tion in General Relativity. In this approach a small astrophysical object follows a trajectory 

consisting of a chained series of bounded geodesies (orbits) in the outer region of a Kerr Black 
&X), Hole, representing the space time created by a bigger object. In our paper we study the entire 

class of orbits, both of constant radius (spherical orbits), as well as non-null eccentricity or- 
bits, showing a number of properties on the physical parameters and trajectories. The main 
l^r^ ' result is the determination of the geometrical locus of all the orbits in the space of physical 

■ parameters in Kerr space-time. This becomes a powerful tool to know if different orbits can be 

connected by a continuous change of their physical parameters. A discussion on the influence 
of different values of the angular momentum of the hole is given. Main results have been 



o 
o 



' obtained by analytical methods. 

PACS numbers: 0420q, 0430D, 0470B 
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H ■ 1 Introduction 



Orbits (geodesies with null or non null eccentricity) in the outer region of Kerr space time play 
an important role in the description of important astrophysical phenomena near black holes. This 
space-time, characterized by the parameters M and a interpreted as the mass and angular momen- 
tum per unit of mass of a rotating hole, can be described in the outer region (limited by the outer 
horizon and the asymptotically flat region) by the Boyer-Linquist coordinate = {t, r, 9, (j)} . 

Geodesies are characterized by four constants of motion, that are E, to, L and C, representing, 
respectively, the energy as measured by an observer at infinity, the mass of a test particle, the 
component of the angular momentum along the axis of symmetry, and the Carter's constant, that 
can be interpreted as the parallel component of the four momentum when the particle crosses the 
equatorial plane of the hole. Sometimes, it will be interesting to use Q instead £, where Q = 
PaPpK""^ > 0, where K is the conformal Killing tensor of the Kerr space-time and = dx°' /dX 
is the four momentum of the test particle (remember that Q ^ C + {L — aE)'^). 



*Also at Laboratori de Fi'sica Matematica, Societat Catalana de Fi'sica, lEC, Barcelona. 
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An extended summary of geodesies in Kerr space time can be found in the well known Chan- 
drasekhar book "The mathematical theory of black holes" (1983) [7]. The algebraic complexity of 
the geodesic equation limited the analytical analysis to particular cases (geodesies in the equato- 
rial plane, in the extreme Kerr geometry a = M, etc.). More recently E.Teo [8j studied spherical 
light rays out of equatorial plane finding an analytical expression for the amount of azimuth for 
a complete latitudinal oscillation of a spherical orbit (the Wilkins effect). This work extends an 
early result of D.C. Wilkins [9] devoted to particle spherical orbits in a = M case. 

More recently, Drasco and Hughes [4j have studied the fundamental orbital frequencies in the 
" r" , " " and " (j)" motion of a test particle using Mino time [5] . This work improves an elegant 
proposal of Schmidt [5], based on the "action angle variables" formalism [6]. 

During the last decade, several works have been devoted to the study of one of the most 
important problems in gravitation, where geodesies play an important role, that is, to compute, 
according to the General Relativity, the gravitational radiation that one can expect from different 
astronomical sources. One kind of these sources are the so called extreme mass ratio inspiral, 
EMRI's: These binary systems arc formed as a result of the capture of a compact stelar remnants 
by supermassive black holes in the nuclei of galaxies. The dynamics of such systems are (almost) 
purely gravitational, accurately modeled as a point particle of mass m, representing the star, 
following a trajectory in the background space-time created by a Kerr Black Hole. If the ratio 
to/A/ is infinitesimal the particle moves along a geodesic. In this system, the particle, having a 
finite mass, radiates gravitational waves. In such a case the trajectory deviates from a geodesic, 
but as the time scale of the orbital evolution is notably smaller than the typical time scale of orbits, 
the particle (approximately) slowly passes from one geodesic to another, conforming the so called 
"adiabatic approximation". Then, radiation reaction effect is characterized by the time evolution 
of the parameters E, L, Q of orbits. 

In order to compute this time evolution of parameters different approaches have been built: 
using of post Newtonian methods, conservation laws, direct computation of self-force and time- 
domain numerical simulations. See for example, among others: Sago and alt. [3] which have ob- 
tained an analytic formulae for the change rates of the energy, angular momentum (using balance 
argument) and Carter constant (based partially on a Mino proposal [2] using the radiative field) (see 
references therein), or S. Drasco and S.A. Hughes [1], that have recently studied the gravitational 
wave snapshots of generic extreme mass ratio inspirals, or the 2007 work "Improved approximated 
inspirals of test-bodies into Kerr black holes" of Gair and Glampedakis [16], etc. . 

In all these papers adiabatic approximation is implemented, and thus, an accurately knowl- 
edge of orbits in Kerr black hole is needed. As an interesting example, we focus on the work of 
different authors that have contributed to the study of the gravitational radiation effect on the 
spherical orbits showing that one spherical orbit remains spherical after radiation in the adiabatic 
approximation (see [TU], etc.). 

In this context, our main goal is to carry out an exhaustive analysis of orbits of test particle in 
the outer region of the Kerr space-time improving our theoretical knowledge, but trying to show 
how one orbit is related to each other, in other words, if they can connect through a continuous 
change of parameters. Wc extend the analytical results on the equatorial plane to the whole 
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outer region finding the set of values that constants of motion can reach and showing the range of 
coordinates that we expect for an allowed set of constants of motions of physical null and non-null 
eccentricity orbits. This analysis allow us to have a global view of all orbits that we show in a 
formal three-dimensional (normalized) parameters space , Q, L, and where 

Q - L - E 
g EE ^ ^ , L= , E=—. (1) 

In this parameter space each orbit is represented by a point. The whole set of orbits is confined in 

^ 1/2 

a compact volume (see figure ^ that partially diverge when E —t 1 and vanishes when E < Za 
(see equation ()62|) for Za,b definitions). Spherical orbits are represented by an upper face (stable) 
and a lower one (unstable). Slices of this volume with E = constant come first {E slightly less 
than 1 ) a pseudo-rectangle (upper side to stable, lower side to unstable, the two bases for orbits 
on the equatorial plane) , and after (za < E < z^^^'^) a pseudo-triangle ( the two sides of 
spherical orbits joints in one vertex when z = zi, see (j62p ) showing that this triangle vanishes when 

— 1/2 

E < Za (see section 4). This representation in the space of parameters allow us to study the 
different paths that a particle can follow going through different geodesies after a continuous change 
of constants of motion. As an specific application, we obtain the necessary and sufficient condition 
that evolving parameters have to accomplish to go from one spherical geodesic to another, showing 
its limits of application. This result can also be found in |10| and |12j in different scenarios and 
using different techniques. 

In section II, III and IV wc analyze the geodesies out of the equatorial plane in Kerr space time 
using well known results in literature I7|,[S],I3- In section V and VI we construct the abstract 
3-space of parameters showing the geometrical locus of spherical and non null eccentricity orbits 
in that space, respectively. Some proofs can be find in Appendices I and II. 

First order differential equations for geodesic motion of massive particles. 

Using the Hamilton- Jacobi method. Carter found the equations of motion for test particles (see 
for example [H]). In the Boyer-Linquist coordinates these equations are 

= ±{[{r^ + a^)E - aL]^ - Alm\'~ + (L - aE)^ + Cy/^ (2) 



where 



p'- = ±{C-[{m'^E')a^ + -^]cos'0y/' (3) 

oA sm 

p^^ = -iaE-^) + &^±^^l^ (4) 

OA sin 9 A 

odi / • 2/, ^^ {r^ +a^)[E{r^ +a^)-aL] . . 

p^— = -a(aEsin^ 9 ~ L) + 4 -. (5) 

dX A 



p2 ^ + c?cos^9, (6) 
A = -2Mr ^ = {r -r+){r -rJ), (7) 



'r± 



M±x/AP ~ a^, < r_ < il/ < r+ < 2M. (8) 



Without loss of generality, in the intermediate computations we are going to work with the 
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following dimensionless constants of motion, often used in the literature (see for example [7]): 

r; = ^, e = (9) 

E E^ 

where £ = C/m'^M'^, considering that z = for photons, and z > for massive particles. In some 
cases, results exhibit a simpler form in terms of 

Q^Q/E^ (10) 

instead of 77. We have to remember that, according with this definition, 

V^Q-i^-af- (11) 

Replacing these three quantities in the equations of motion and taking into account that for 
massive particles A = r/m, where r is the proper time, the equations become 

= ±z-iVR (12) 



^ ±2z-^Ve (13) 
ar 

- 2z-2[(- )^+ ] (14) 



df ^'1-^2 ^ 

J dt 



= z--{a[i^a{\- + ^ ^}. (15) 



where 



n A 

R = -x(x-l)(f+ 7)'-y{^ -[z^{z~l)x]} (16) 

a; — 1 4 X — \ 
e = (l-M')?y-M'e'-A^'a2(2-l)(l-^2)^ (^7) 

p2 = 4x2 + aV^ (j^g) 

A = Ax^ -Ax + (19) 

r - i _ T 

a; = , u = cos0, t= , t = . (20) 

2Af ^ ' 2M' 2Af ^ ^ 

The allowed values of the independent constants of motion ^, r\ and z are only limited by 
conditions i? > and O > 0. Equations ([T^ and pS]) give us the variation of </> and t on proper 
time, once we have chosen allowed values of these constants. 

From now on we will drop the bar symbol from a, A and p considering that we will ever work 
with non-dimensional quantities. 

As it is well known, these equations of motion are valid for all values of a. However, in our 
work we have differentiated a = case (The Schwarzschild space-time) from the < a < 1 case. 

Furthermore, we are interested in the bound geodesies (z > 1) in the outer space-time of the 
hole, that is, a;+ < a; < 00. 



2 Theta-motion for bound orbits 



It is well known that, without loss of generality, we can study all geodesies in Schwarzschild (a = 0) 
space-time in the equatorial plane (Q = 7r/2, dO/dr = 0, 77 = 0). 
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In order to analyze the whole set of geodesies in Kerr space-time < a < 0, the 6'-niotion 
equation sets important restrictions on the ^ and rj values. 

Considering z > 1 case (we include z = 1 as a special case) condition > holds if and only 

if [2 

?7 > 0, (21) 

-M- < ^^<^^-, (22) 

in accordance with (|17p . where 

= + + Vie +V + a"? - 4a27?} (23) 

za^ 

< fi^_< 1, 
= (z — l)a^, 

In the z = 1 case we find that if r/ = 0, then /i = (^^ = n/2) or ^ = and < /^^ < 1. 
Moreover, if > geodesies may be found. 

In the z > 1 case, r/ = ^ M^(C^ + a^(l - /i^)) = 0. It implies: /i^ = ^ = |, or 
^^-|-a^(l — /i^)=0^^ = and 9 = 0,tt. Besides, if 77 > geodesies exist. 

According to (|23p . 1 = 1, if and only if ^ = and r] > a^. It means that only the trajectories 
with = can reach 9 ^ Q,tt, that is, the symmetry axis. This result can be found not only in 
particle geodesies case but also in light rays geodesies (see [5]). 



3 Locally Nonrotating Frame 

There's a condition that arises when we consider that the photon's or particle's energy measured 
by any inertial observer cannot be negative or zero. We know that if one inertial frame measures 
positive energy then all inertial frames must measure positive energy. We will use an observer who, 
in some sense, rotates with the geometry, the so called 'locally nonrotating frames LNRF" (see 
[13j for definitions). Thus the observer world line is r = constant, 9 ~ constant, (p = ujt + constant 
where w — —gct,t/gtj)phi- The orthonormal tetrad carried by this observer is 

e. - (J_^i/2d_ 2Mai_d_ _A.i/2d_ .l^l/2 5_ _ ( ^l/2 d 

' " at ^ (Ap2A)i/2a0' ""'-^^p^' dr' '~^p^' 89' ^"^Asin^fl^ df 

A = {r^ +a^f -a^Asui^9. (25) 

In this frame the measured energy is 

Elnrf = -P"et^ ^B^{E- MwL), (26) 

where J5 is a non-diverging function and w is 

yj = ^ (27) 

9vv A/[(4x2-t-a2)2-a2Asin2 6']' ^ ' 
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It will be interesting to find bounds on jj^ . These are 



Mwm ~ Mw ~ MWM 

where 

Aax 



^ > t:^ > > 0. (28) 



Wm 



M(4a;2 + a^)^ 
a 



It is important to remark that 



wm{x+) = Wm{x+) = — — , „, > 0. (29) 

zM (1 + V 1 — a ) 



The requirement Elnrf > is equivalent to £■ — MwL > 0. Therefore 
E > ^ ^ < < ^ 



Mw ~ MWm 

E < ^ ^> — — > — ^ > ^ L < 0. (30) 
~ Mw ~ Mwm 

In order to apply this rule out of the equatorial plane, (?; 7^ 0), we use the relation (pS)) showing 
that takes values bounded by two 6'-non-depending limits. 

In the Schwarzschild case = and then Elnrf > i? > 0. 



4 Spherical orbits 



Orbits with constant Boyer-Linquist coordinate "r" {x = constant) require that i? = 0, and 
R' = 0, where R' = dR/dx. In the literature one refers to these orbits as 'spherical' or 'circular' 
orbits indistinctly. The analysis of spherical geodesies stablishes a basis to classify the whole set 
of geodesies. We focus our work on the outer region of Kerr space-time x+ < x, as we have said 
before. 

Solving simultaneously R = and R' = equations we obtain general expressions for and 
r]s , in terms of Xg , where a and z are taken as parameters (the subscript s means that these values 
correspond to spherical orbits) . The solutions for are 

_ 4xs'^ - ± Asy2zr^2(l - z)xs + z 



(2x, - l)a 



(31) 



while the corresponding solutions for rjs are 

4x2 



'7sl,2 



{2xs - 1) a2 
4xs(4a;^ - Sx^ 



z{2xs - 1)(8< - mxi + 8x, - a^) 



5x^ 



T 2A,\/2^V2(1 - z)xs + z 



(32) 



If 1 < z, the necessary and sufficient condition for these solutions to exist is: 
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Note that Xg^^^ is an upper bound to the radms of spherical orbit. This hmit is independent of a, 
as long as a 7^ 0. 

From now on. we will focus on the following intervals for x, a and z: x+ < x < x^^^^, 1 < z, 
< a < 1 that we will generically call T), the domain of variation of these quantities. 

In order to ensure the existence of V it is necessary that x+ < Xs^^^ ■ This implies an upper 
bound of z, that is 

1<^<2:._, (34) 

where 

Zg_^^l + —L=, (35) 
V 1 — a 

(a = 0, Zs^^^^ ~ 2; a ~ 0.8, 2i5,„„^ = 2.66..; a — > 1, Zs^^^^ oo). 

If 1 < z, the 0-motion analysis has shown that only non negative values of rj need to be 
considered. Then circular orbits occur when 7751,2 > 0. 

After a tedious calculus we can see that 

risi{xs,a,z) <Q (36) 

in T) (see Appendix I). Therefore spherical orbits in Kerr space-time (0 < a < 1) can be foimd if 
and only if i]s2 > 0. 

A. The {ris2,is2\ solution. 

Firstly, we apply the results of section 3 to this solution. We can check that 

- jt;-, - M^' ^^^^ 

(if X = .T+ we get the equality). Therefore, according to ([30|): No spherical orbits exist if £^ < 
and all spherical orbits with E' > have Elnrf > 0. 

The analytical complexity of these functions doesn't allow us to have a complete knowledge of 
them in a direct way. While 7752 has three extrema in T> (every extremum is a maximum, minimum 
or inflection point depending on the a and z values), 2^2 1 defined in (jlip . 

Qs2=Vs2 + {a-L2? (38) 

only has one and shows other interesting properties that we will discuss below. 
The function R can be written as follows 

R = A{S,+Bf + C{Q + g), (39) 

where A,B,C,g are 

A^—, B{x,a)= , C(x,a) = - — , g{x,z)^izx^. (40) 

Thus R must be considered as a function of x,a,C^, i.e. R{x,a,C^), where C" = {^, Q, 2}. 
According to this we define 
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For each set of allowed values of C" we can find a solution x = x{t, C*, a) solving equation p^ . 
Thus R = R{x{t, C\a),C^,a). 

We arc now interested in a family of solutions such that the parameters vary according to a 
specific law, defined by the functions C"(s). assuming that for each value of s we obtain a geodesic 
[T2] . The solution of (fT^ now takes the form a; = x{t, C"(s), a), and thus 

R^ R{x{T,C'{s),a),C^{s),a). (43) 

If we produce a small variation (5s then the quantities R and R' vary according to 

_ . s ||c-. (45) 

where x = dx/ds and = dC^/ds. 

We apply this results in order to investigate the properties of the first derivatives of £,s2ixs, a, z) 
and Qs2ixs,a, z). 

Spherical orbits imply a; = s = Xg. This means x = 1, thus, ^(s = Xs) = S,s2{xs, z,a)^ 
Q{s = Xg) ~ Qs2{xs, z, a) = r]s2{xs, z, a) + (a — ^32(2;^, z, a))^ and z — 0. For every allowed value of 
s = Xs we have spherical orbits, therefore R{xs, ^^2, Qs2, z, a) = R'(xs, ^^2, Qs2, -z, a) = 0. As we go 
from one spherical orbit to another varying s, hence R{xs, ^s2, VS2, z, a) = R'{xs, Cs2, Vs2i z, a) ~ 0. 
Using these results on (jS]), (|44|) and ([45|) we obtain 

= 2AiU + Bs)B'^ + C'^{Q,2+gs) + C,g'^, (46) 

= 2^(e.2 + Bs)e:.2 + C^.Q;2 (47) 

= R': + 2ABU's2 + 

+ C'M2, (48) 

where for example = B{x — Xs), and ^^2 = d^s2/dxs = Cs2- ™ accordance with our previous 
notation. 

Theorem I: £^g{xeia, z), a, z) = if and only if Q'g(xe{a, z), a, z) = 0, where Xs ~ Xe{z,a) is 
a solution of ^'^{xs, a, z) — Q'g{xe{a, z), a, z) = 0. In other words, there where ^si^s, a, z) has an 
extremum, the function Qs{xs, a, z) must have an extremum, and vice versa. 

Proof: Considering that ^^2 and Qs2 satisfy i?'(^s2, ?7s2, a^s, a, z) = 0, then, according to (|T7)) . if 
£.s2 — then = since Cg = — As/4 < (a;+ < Xg < Xg^^^)- Conversely, if = it implies 
that " since the other possibility ^^2 = —Bg, g'g ~ and Qs2 = — ffs (to ensure R = R' = 0)is 
impossible due to > and Qs2 must be greater than zero. 

In Appendix III we have proved that this extremum exists and is unique. We call it Xe{a,z) 
but we are not able to found his analytical expression. 

From equation ^ '^c can obtain 

8a;„ 



{32x1 - 36x2 i2x^ _ a2)2 



[2xe{2xe - l)(32a;^ - 44a;2 + 20xe ~ 30^) 



3 



AiV2xe]. (49) 
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Using this expression we find 

^e2{Xe,a) = ^s2{Xs = Xe,a,Z = Ze) (50) 
Qe2iXe,a) = Qs2iXs ^ Xe,a, Z = Ze). (51) 

where the functions ^e2ixe,a), Ze{xeTa) and Qe2ixe,o-) represent the values of minimum of ^^2, 
maximum of Qs2 and Zg when the extrcmum is in x^. 

Theorem II: The fimction R"{^s2, Qs2, Xs, a, z) vanishes if and only if t ^32(^3, a, z) ~ Qs2ixs, a, z) 

0. 

Proof: The right to left implication is obvious using (|48|) . The left to right implication can be 
proved considering that the determinant of unknowns ^^2 ^.nd 2^2 the homogeneous system of 
two equations ((47l) and (|48p . considering that now R" = 0, is 

A, 



Det = 2As{U + B,)C', - 2A,B'fi, = -f^ -2x,{2x,z ~ 2xs ~ z). (52) 

It is clear that Det only vanishes 'mx ^ x+ and x ~ Xg^^^ ■ Then the only solution is Cg2ixs, a, z) = 
Q'g2ixs,a, z) = in x+ < x^ < Xs,„^^ ■ 

As it is well known stable (unstable) spherical orbits occur when R"{x, a, Qs2,£,s2, z) < 0(> 0). 

Theorem II implies that R" = only where Q'g2ixs, a, z) — £^'g2{xs, a, z) — 0, that is. where 
X = Xe{a,z). It is not difficult to prove that spherical orbits arc stable when x^ > Xe{a,z). 
Conversely, if Xs < Xe they are unstable. 

Theorem III:The extremum of ^^2 and Qs2 are maximum (minimum) and minimum (maxi- 
mum) respectively in T). 

From ([47|) . and considering the point where ^^2 = Q's2 = we obtain 

W2 ^ '2As{^g2 + Bg) < 
Numerical calculations show that this point is always a maximum of Qs2 and a minimum of ^^2 • 
The values of these two functions in xj^ and Xs^^^ are 



Qs2{x+,a,z) = -(1 + Vl - a?)z < 0, 
Qs2{xs^^,,a,z) = =-j—j[z{l-a^)+a^]<0, 

U{x+,a,z) = ^ii±^^iHZ) > 0, (54) 

a 

U{xs^^^,a,z) = [z(l-a^) + Q^]+a>0. (55) 

a[z — 1) 

B U = ^2 = 0. 

It will be important to find the functions Zc{a) and Xec{a) such that if z > Zc then S,s2ixs-, a, z) > 
0, where Xec is the point where = Cs2 = 0, for any value of a. To do this we proceed as follows: 
First, the system of equations ^^2 = Cs2 = is equivalent to 

1/2 



2x 



1/2 



+ 1 + ^ {Aix^^y/^ - A(XeJ-l/3 + V2) ' 



(56) 



1 (4 x,J + g^) (8 Xec^ - 12 Xec + 2 g^xesc + g^) ^g^^ 
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where 

A{xec) = ll\/2 - 27V2xec + 3^/27 - 132x^0 + 162x2^. (58) 

If we use ([55)1 in ([57)) wc obtain a function Zc{xec) that combined with is a parametric 
expression of 2c(i) (see figure (Il]),a), that is, for every a we find Zc. Then, using (|58[) (see figure 
([1]), b) we can find Xec value. 




Figure 1: Plot of Zc{a) and Zc{xec) 

Numerical examples can be done: a = 0,Xec = 3,Zc = 1.125; a = 0.8, Xec = 2.784, Zc = 1.133 
and a 1, Xec = 2.637255281, Zc = 1.138373645. 

As an interesting application of this result, wc can prove that d(t)/dT > for all spherical orbits 
provided z > Zc- According to p4|) 

z'^jj^dd) , 1 . ^ Aax , , , 4aa; , , 

-^Tf-^T^~ A + ^ ^ - A + 

But 



Therefore 



„2 a l + v2\/— (— 2a:: + 2za; — z)a; 

z^p^dd) , 1 4aa; , 

^df = (T^-A)^- + ^>^- 

Thus, we can conclude that: 

Theorem IV. For every value of a, there exists Zc, defined by ((57|) and (|56p . such that if z > Zc 
then > in I? and thus d<t)/dT > 0, i.e. the spherical orbits co rotate with the hole. 

C Qs2 - (a - 62)' = Q's2 - [{a - U?]' - 0. 

In order to improve our knowledge of Xg, we can find where both Qs2 ~ [o- ^ ^^2)^ and ~ 
[(fl ~ ^s2)^]' 0, or alternatively ?ys2 ~ Vs2' = 0. These conditions arc equivalent to 
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Xea,b = ^{3 + Z2T[{^-Zi)(3 + Zi+2Z2)]h- (63) 

where 

Zi = l + (l-a2)*[(l + a)* + (l-a)*], 

Z2 = {3a^ + Zi^)i. 

Functions Xgj, f, where first presented in |13j (see equation (2.21)in that paper). We now show some 
values of these variables: 

• When a = 1 then = 3 and Zf, ~ 1.080 while Xea = 1/2 and Xeb = 9/2, 

• When a = 0.8 then Za = 1.298 and Zb = 1.086 while Xea = 1.453 and Xeb = 4.216, 

• When a ^ then Za, Zb 9/8 while Xea, Xeb — > 3. 

Moreover, Xea,b give us the range of variations of Xe, that is Xebia) < Xe{z,a) < Xea{o) obtaining 
non-negative values of r/s2. 

Finally, we can verify that 1 < Zb < Zq, see fig ([2]). 

1 ,13 
1,1275 
z 1,125 
1,1225 



0,0 0,25 0,5 0,75 1,0 

a 




Figure 2: Dependence on a of z;,, and Za- 
The behavior of these functions near a = 0+ is 

2a = 8^32""^^^- 
this implies that Zb < Zc < Za (the equality holds in the limit a — > 0). 
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D Classification. 

Until now we have studied the functions 2,2 and ^,2. We conclude this subsection analyzing 
under which conditions Qs2 > (a~Cs2)^ (<;=^ 7752 > 0) holds. We call x^, i, j, .. = {1, 2, ...}, Xj < Xi+i, 
the points where Qs2 ~ {a — Cs2)^- We will depict these two functions Qs2 and (a — ^s2)^ in four 
possible scenarios (see figure 

• 1)1 < z < Zf,. Condition only holds in two disconnected intervals, a;-|_ < xi < x < X2 < x,, for 
the unstable spherical orbits (in which Qs2 is an increasing function) and < Xe < X3 < x < 
Xi < Xsmai that corresponds to the stable ones (in which Qg2 is a decreasing function), 

• II)z = Zf). The two parts are just connected in Xg, that is X2 = Xe ~ x^ of the previous case, 
and we can continuously go from the stable to the non-stable orbits, 

• III) zi, < z < Za- The condition holds in the whole interval x+ < Xi < a: < X4 < Xs^^^ 
containing stable and unstable orbits (now X2 and X3 doesn't exist). 

• IV) z = Za We only have one spherical geodesic in Xe- 
In figured! we can see ^^2 in these cases, including z = Zc- 




Figure 3: Plot of Qs2 (continuous line) and (0—^32)^ (dashed line) corresponding to /; //; ///; IV 
cases 

E Equatorial Plane. 

As it is well known, a necessary condition to have orbits lying entirely in the equatorial plane 
is 7] = Q — (^ — aY = 0. For spherical orbits the condition becomes i]s2 ~ 0. 
Conditions R = R' ^ are now (77 = 0) equivalent to 

[16 x,^ - 56 x/ + 64x^3 _ 6 + 4) x^^ + 5a2^^ ^ aAs^/257] x^ 



16 Xs'^ + 8xs — a?Y 



^ ^ _ -q(12x^^ - 8xs + a^) T 2 AsXs72x7 ^^^^ 

8 Xs"^ — 16 Xs^ + 8 Xs — 
Qs± - (6±-a)'. (65) 

These equations arc equivalents to (2.12) and (2.13) of [T3] . 

Each of Zs± solutions has a maximum in the interval 1 < z < 00, that corresponds with 
z^M = Za,b, x±M = Xea,b defined in (|63p . Beyond these points they decay asymptotically to one, 
as we can see in figure ([5]). 

As you can see in the last figure, if 
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Figure 5: Graphic z(xs), rj = 0, a = 0.8 

• 1 < z < Zh we have four spherical orbits. 

• z = we have three, one of them in Xeb such that {zb,Xeb,<i} are solutions of (|62p . ((63)) 
equations {zi, < Za), 

• Zb < z < Za we have two, 

• z — Za 'wc have the last one, such that {zq, Xea-, a} are its solutions. 

A specially simple expression of the function Zg± (xg , a) can be obtained for z = 1 , and z = 9/8 
values, that are well known in the literature. 
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a{z = 1) = J4{xs - 2^2^+ 2)a;, = 2Jxs{^/^s- V2)'^. (66) 



a(z 9/8) = ^^^y^-7a;,2 + 18a;, + 81 ± 4(9 - Xs)Vx{Q - 2a;,). (67) 



9 

(See [H] and [7] for details). 

F Schwarzschild circular orbits. 
Conditions R — R' — imply 



= ±^ (68) 

Xs - 1 

= 2(a;, - ir ■ ^^^^ 

Q.± = e±. (70) 

Functions and have a maximum, ^s-^/ = ^^^^j -^sa/ = 9/8 respectively, in = 3, and 
a minimum, ^s+m = 

Stable circular orbits can be found when 3 < a;,, < cxd, and unstable 2 < x < 3, where 
Xs = 2, 2, = l,^s± = ±4 corresponds to the inner unstable circular orbit. (See [13] and [7] for 
details). 



5 Spherical orbits in the Q — L — E space 

In this section we are interested in representing "outer" spherical orbits in a Q, L and E ([T]) formal 
three-space, ( Q ~ z~^Q, L = z^^/^^ and E = z^^l"^)^ analyzing their properties, as we have 
mentioned before. This task can be achieved due to the special shape of the geodesic equation 
that, at the same time, derives from the existence of two killing vector fields and one conformal 
killing 2-tensor. As we have seen t,si and Qs2 depend exclusively on x = r/2m coordinate. This 
fact allow us to construct a formal three space in such a way that the whole set of spherical orbits 
will be represented by a 2-surface, say E, in this space. In the next chapter we will see that orbits 
with non null eccentricity can be represented in the same formal space. Therefore, we will achieve 
a global and careful (with accurate analytical control) representation of the whole set of orbits in 
Kerr space-time. 

A. Each point of this space is characterized by their coordinates = {Q, i, E}. 
The geometrical locus of the spherical orbits in this space is a 2-surface S defined by: 

Q = z^^Qs2ixs,z,a), 
L = z^i^s2ixs,z,a), 

E = (71) 



limited by 



Q = {L-aEf, 



(72) 
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that is, the surface representing the equatorial plane in this space (77 = 0), and < E < 1. Finally 
Xs and z (z = E~'^) have to be considered the two parameters constrained by 

x+{a) < Xs < Xs^^^ (z), 

I < Z < Za{a). 

where Za is defined in ([62)) . Wc will look at this representation a(0<a<l)asa constant (for 
example a = 0.8), discussing its variation in the results and including z = 1 as a limiting case. 

Moreover, this representation is not sensitive to initial conditions. Then, one point represents 
spherical orbits but, in some circumstances, the same point represents the associated geodesies 
(see next section for definition). 

It is well known (see [14) for example) that singular points occur when two minors of the matrix 

dQ dL dE _ 
dxs dxs dxs 
dQ dL dE 
dz dz dz 

vanishes. This condition holds if and only if = Q's2 = 0- The solution of these equations is 
Xs = Xe{z, a). Replacing this solution in ()7ip we obtain the set of singular points of S. 
B. 

Spherical orbits can be represented using a parametric plot of equations ([7T)) . The range of 
variation of parameters Xs and z are such that the surface S is limited by £' = 1 plane and the 
surface Q = {L — aE)^ as is shown above. In figure © we depict the intersections of these surfaces, 
including slices obtained cutting S by planes of z = constant, as we explain below, for a particular 
values of a. 

In this figure we can sec how S is in fact a surface folded by the line from point A to point 
C3 made of singular points. The equation of this line has been found in ([49[) . ()50)) and ()5ip . Each 
point represents the spherical orbit with the maximum value of Q and the minimum of L for a 
fixed value of z and a. 

Unstable spherical orbits are represented by the points of the lower pseudo-trapezoid face with 
curved sides (that we represent by a widehat over the extrema of the piece of curve). These sides 
are: 

• A, El and Cz,E^2, unstable circular orbits on the equatorial plane, 

• A, C3, singular points of S, 

• El, i?32, unstable spherical orbits with E = 1. 

Stable spherical orbits form a partially unlimited face with sides: 

• A, C3, singular points of S, 

• A, E2 and C3, E^i, stable circular orbits on the equatorial plane. 
The apexes of these partially non-limited volume are: 

• A the circular orbit in the equatorial plane with smallest value of E for a fixed value of a. 
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• El and E32 circular unstable orbits in the equatorial plane with E = 1, 

• C3 circular orbit in the equatorial plane with smallest value of E and L < 0. 

• E2 and E31 limiting points representing stable circular orbits in the equatorial plane when 
E 1 

We complement our knowledge of S cutting slices of E ~ constant{<^ z = constant). Then, 
the spherical orbits could be represented in each slice in a parametric way by ((7T|) where z takes a 
constant value. We include the intersection of each slice with the surface of the equatorial plane. 

The complexity of the behavior of these parametric functions forces the study of their particular 
shape for the cases studied in table 1. Relations listed in the previous chapter allow us to better 
understand those particular shapes. 




Figure 6: The geometrical locus of spherical orbits in Kerr space time. Slices are made for I, z = 
1.05; //, z = Zb = 1.086; ///, z = z^ = 1.133; IV, z = Za = 1.298 

• Case 0: z ~ 1. Ei, E32 unstable spherical orbits with E = 1 

• Case I: 1 < z < Zb.Di,D32 unstable spherical orbits, D2,D^i stable spherical orbits, 

• Case II: z = z^.Ci, C3 unstable spherical orbits, C2, C3 stable spherical orbits E ^ Zj^ 

• Case III: z^ < z = Zc < Za Bi, unstable spherical orbits, B2,B3 stable spherical orbits 

E = z-'l\ 
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Figure 7: L — E plot of circular orbits on the equatorial plane, Q = {L — a)^. 
• Case IV : z = Za Last stable spherical orbit. 

In order to complement the previous plots, we add the projection on the E,L plane of the 
spherical orbits on the equatorial plane. We represent these orbits on a parametric plot using (|64p 
into ([71]) . See figure ([7]) plotted for a = 0.8, that arc the projections of A,Ei, A,E2, C3,£'3i and 

^3 , E32 ■ 

As a decreases both peaks of this L — E plot become more and more symmetric until we arrive 
to a = 0, where they become completely symmetric (as could be expected from the spherical 
symmetry of Schwarzschild black holes) 

C As we have seen in section D of the previous chapter, in the general case I, the function Qs2 
has four values of Xg, say Xi i £ {1, 2, 3, 4}, such that Qs2 = (a ~ ^s2)^, for a fixed value of a and 
z. The plots have been obtained in all cases increasing the value of the Xs parameter, first from xi 
to X2 (in figure [6] from Di to -D32) for the unstable branch, and from X3 to X4 for the stable one 
(from Z?3i to 02). In the other cases the end of the unstable and the start of the stable branches 
is Xe which traces the singular points of the global surface. 

We can now show the relative values of £,s2{,Xi, z, a), that is 

Xi £,s2{xi,z,a) = ^i 

I Xi < X2 < Xs, < Xi 6 < 6 < < 6 < ^4 
II Xi < X2 = Xz < Xi 6 = 6 < < Cl < ^4 
/// Xi <Xi < ^1 < ^4 

IV xi = Xi < = ^4- 

Moreover, since is a continuous function, it has a minimum in x^ and takes positive values 
in the extrema of T>, x^ and Xs^^^ , from the above table, we can infer that ^s2 = in two 
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points xi {S.s2i^h '^jO.) < 0) and Xm {£,'s2{xm, z, a) > and xi < Xm) if 1 < 2 < Zc- Thus: i) 
Xi < xi < X2 < X'i < Xm < xa in the case I and II, ii) if < z < z^., then x^ < xi < xi and 
Xe < Xjn < X4, iii) if z = Zc then — xi — x„i- 

D Continuously moving in a geometrical locus of spherical orbits. Let us choose P, a point of 
E. The holonomic base of the tangent space in P associated to curves on S of = constant and 
z = constant is 

X, = A =Q,A+L,A+44. (74) 
9z dQ dL dE 

oa;s dL 

where Q' and means partial differentiation of Qs2 with respect Xs and z respectively and 
S:,^,^ is the holonomic base on P of Q,L,E space. 

oQ dL dE I 7 

As a consequence of gravitational radiation, in the adiabatic approximation, a test particle can 
follow a trajectory that can be imagined as a chain of orbits in such a way that the particle remains 
several periods on each orbit before changing to another orbit. We define t as an average over 
several orbital periods. According to this, the constants of motion {Q,L,E) associated to different 
orbits vary according to a law that evolves dQ/dt,dL/dt,dE/dt quantities. In fact, the solution 
of the first order ordinary equations is a curve in our formal space such that 

V = ^A + 4A + ^A (76) 
dt dQ dt dL dt dE' 

where V = d/dt, is its tangent vector (' = d/dt). The necessary and sufficient condition to ensure 
us that a particle in spherical motion remains in a spherical orbit under the influence of radiation 
reaction is 

e^.kV'X^jXj' = ^ (77) 
^ dQ QjJ-Q^dE ^,dL 
" 1^ d^ + ^dt' 

where e^fc is the Levi-Civita tensor. Hence, equation ([75]) becomes [H], (see also equation 3.25 of 

m) 

dQ ^ ^2x. + ^^2xA-z-2xs + 2zx,) dE _ J_L 

dt (2xs-l)Vi dt dt^ ^ ' 



If the motion lies entirely on the equatorial plane then equation (|78l) , taking into account ((72|) , 
becomes 

dE J di 



The result is 



—L^ = —E^. (80) 
dt dt ^ ' 

dL dE 1 , . . 2^/^x^A , , , 

— = -7-: —Ux^+a^+ , )■ 81 

dt dt(2x-l)a^ ^-a;(2za;-2a;-z)^ ^ ' 

Using basic instruments of classical differential geometry we can find the necessary and sufficient 
condition for a stable spherical geodesic be able to achieve a non null eccentricity orbit following 
V direction, that is, eijkV^XxJ Xz'' < 0. This result apphes if the start point is an unstable 
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spherical orbit and we want to achieve a non null eccentricity orbit in V direction. In general this 
construction allow us if whether or not a specific curve, with V as a tangent vector in each point, 
is inside this partially unlimited volume of orbits in Kerr space time, and therefore if this is able 
to describe the path of a radiating particle in the adiabatical approximation. 



6 Non-null eccentricity Orbits in the Q — L — E space 

In the previous section we have constructed a space where we have placed the spherical orbits, 
now we are going to study all orbits and we will see that spherical ones delimitate where the rest 
of the orbits are. 

A The U potential. In order to analyze the non null eccentricity orbits in relation to the 
spherical ones we could use the potential technique factorizing R as follows, 

R = T{U-Q), (82) 

where T and U are defined below. The complexity of U advise us to use 77 to factorize R. It will 
be easy to translate the results in terms of Q and L as we show below. 
According to ((39l) we have 

R = j[U-v], (83) 

u = — -[z-iz-i)x] — — + 7'- 

X — 1 A X ~ I 

4a; {4(z - l)a:3 - 4zx2 + [^2 _^ (z - 1)0^] X - (C - a)2} 

Therefore, considering that T{tj — Q) ~ ^{U — rj) , wc obtain 

U = zU-z{C-af, (85) 

T = ^, 
4z' 

Q = (86) 
Orbits exist only if i? > 0. This implies 

U >r]. (87) 
For fixed values of ^ and z parameters, and a, U depends on x. 

• X x+ then U +00 

• a; — > 00 then U —00 

• U is continuous in the domain: 0<a<l; 1 < z, < x < 00 

• Using Descartes's rule on the third-degree polynomial part (see(|84l)). one can see that, for 
the values of a, z and ^ that we are interested in, there are always 3 changes of sign. Thus 
there's the possibility of having at most 3 real roots, and wc know that there will be at least 
1 root, because this function, being continuous, goes from 00 to —00 inside that region. 
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Figure 8: Relationship between Q — L plot and U potential in two scenarios: at left Lo = 2 
intersects the stable and unstable curves and so we have both extrema of U above the minimum 
value of Q, i.e 1.48 i] = Q). At right Lq = 4 only intersects the stable orbits curve and then we 
only have the maximum of Q above the minimum allowed value of Q, 10.36. Regions where orbits 
exist are filled with vertical lines. 



Extrema of U can't be found analytically. As it has been explained, the conditions for spherical 
orbits are R = and R' = 0. In U potential terms, they are equivalent to: 



7? = U{x,a,^,z) (88) 
U'{x,a,^,z) = (89) 

The set of values of ^, 77, z, x, a which correspond to spherical orbits arc related through the func- 
tions and 7?si,2- We are going to study the extrema of U through these functions proceeding 
as follows: First, one value of oq and zq is chosen. Then, to completely define the U potential, we 
have to choose a particular value of ^, say ^o- 

The two functions S,si,2 coincide at the extrema of the interval of a;+ < x < Xs^^^. They 
are continuous functions and each one has one extremum in the interval, a minimum for ^^2 and 
a maximum for ^^i- So they form a oval closed shape. Hence, we can find different scenarios 
depending on 

• ^0 is such that > Csi or ^0 < Cs2- Thus U doesn't have any extremum. 

• ^0 is tangent to the oval shape at ^si maximum or ^^2 minimum. In this case U potential 
has an infiection point. 

• ^0 intersects twice the shape. This implies that two extrema exist in D. In this case, in- 
tersections with curves implies that the corresponding extrema are below the horizontal 
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axis since 77^1 < 0. But for intersections with £_s2 the extremum will be over or below the 
horizontal axis depending on the sign of 77^2 • 

As we have seen, [/ is a continuous function that takes the values +00 and —00^ at and 
a; ^ 00, respectively. Then if it has two extrema, it must first have a minimum, and after a 
maximum, with the maximum being above the minimum. 

In figure ([8]) we can see the more general example where the two possible scenarios are plotted in 
terms of Q and L by using ((85|) , and ((86|) (the form of U is like U except a shift and a rescaling) . We 
have complemented the information comparing these two plots with the corresponding Q versus 
i, that corresponds to the Case 1 (slice _Di, Z52, -D31, £'32 of figure © of the previous section), 
establishing graphical analogies. 

Now that we know the shape of U we can find where the non-spherical orbits are. An orbit 
must exist between two turning points, if we choose a value Qo we see that it can intersect up to 
three times the U potential defining three turning points Xt < Xa < Xp. We can find orbits only 
when Qq value is between both [/ extrema (considering Qo > (io ~ o-Eq)"^ , where Lq and Eq are 
fixed values defining U) and Xa < x < Xp, such that the apoaster is in Xa, and the periaster is in 
Xp. This situation corresponds to the region between the two curves at the Q versus L plots and 
the curves ?/ = <^ Q = (ij — aE)^. Geodesies falling into the hole exist in x^^ < x < xt region 
(associated geodesies |T|), assuming the same values of the physical parameters Qq, Lq and Eq of 
the associated orbits. 

Figure ^ tell us that, starting from a value of Qo corresponding to an orbit, if we increase 
continuously this value we can achieve at most the value Umax in order to have physical orbits. 
But if we decrease this value in the left case of ([8]) below Umin wc will get open geodesies that will 
be swallowed by the hole, but at the right case of \i Qq < {L — aE)'^ no physical orbits exist. 
Therefore, the surfaces representing stable spherical orbits and the Qq ^ {L — aE)"^ surface in 
figure ^ are impenetrable, i.e. can be achieved but not crossed following a continuous change of 
the physical parameters from an interior point of the volume of orbits, while the surface of unstable 
spherical orbits can bidirectionally be crossed. 

B The region where orbits exist. 

If we focus on 1 < X region we can obtain interesting results due to the fact that, in that region, 

<V <U{tz,x,a) < h{x,z) (90) 

since 

U = /.-^^(e + ^)^ (91) 
A X ~ 1 

h = .^[z - (z ~ l)x]. (92) 

X — 1 

If a; ^ 1+ then h +00, if x —^ +00 then h —00, besides if 1 < z < 9/8 has a minimum in 
Xei and a maximum in Xe2, x^i < x^i and h(xei, z) < h(xe2, z) > 0, and if 9/8 < z, thus h{x, z) is 
a continuously decreasing function. In addition to this h{x^ z)>Q-i=>\<x< xumit, where 



Xlimit = z/{z - 1). 



(93) 
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Considering that h{x, z) < in the x+ < x < 1 region, we can assert that the allowed values of x 
coordinate of orbits are bounded by a;_|_ < a; < xumiti^)- The upper limit only depends on mass 
energy ratio. 

C Orbits CO rotating with the hole 

According to figure [SI physical orbits with z > Zc arc characterized by L > {■i^ ^ > 0)- In 
fact, there exists a set of orbits with the same value of L. These orbits are placed in this figure 
in the slice oi L = constant > 0, including a number of stable and unstable orbits. Therefore 
Theorem IV must be extended beyond spherical orbits, concluding that all orbits with z > Zc co 
rotate with the hole. That is, following ([T^ we can check that dip/dr > for all orbits with z > Zc- 

2 df 4-//2 A^^+ A A^^^ A- ^ > 

But, we can always find a spherical orbit such that ^ = ^s2- Thus 



Aax ^ , a . ^ Aax 



a(^l + \/2^- (-2 x + 2zx- z)x^ 



Therefore 



z^p'^dch , 1 c? . Aax , 



7 Summary of main results 

The main goal of this paper has been to provide different tools in order to have a very important 
knowledge of both (unstable and stable) spherical orbits, and non null eccentricity orbits in the 
outer Kerr space time, studying how their physical parameters arc related. This is due to the 
special role played by Kerr orbits to explain astrophysical effects in strong gravitational fields such 
as gravitational wave radiation in the adiabatic approximation. 

To achieve this, we have developed an analytical study of orbits out of the equatorial plane, 
recovering, as expected, classical results of orbits lying entirely on the plane of symmetry including 
the Schwarzschild results and the limit of Kerr metric (a = 1). A numerical approach to obtain 
specifical results about orbits or radiative models can be implemented based on our analytical 
approach. 

1) We have classified spherical orbits in 4 different classes according to the different behavior 
of the physical parameters restricted to spherical orbits Qs2, Ls2 and E. The infiucncc of a in the 
results has been studied, as well as maximum values, domain of existence and correlations of Qs2 
with La2, showing that these two functions have a common extremum in e ( a minimum of Ls2 
and a maximum of Qs2) that, at the same time, is the threshold that divides stable and unstable 
spherical orbits (see Theorem II). 

2) Moreover, using the properties of the extremum e mentioned above, we have found that, for 
z > Zc(a) Ls2 > 0), all orbits (including the spherical ones) have a positive value of L and 
therefore must co-rotate with the hole. 

3) According to this classification we have constructed the physical parameters space (Q, i, E) 
showing that the whole set of spherical and non-null eccentricity orbits can be represented there. 
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We have been able to carry this out thanks to the fact that the Kerr metric has two kiUing 
vectors and one kilhng 2-tensor and the choice of Boyer-Linquist coordinates. In this space the 
geometrical locus of spherical orbits is a 2-surface that looks like two sheets (one for stable and 
the other for unstable spherical orbits), partially matched through a common curve. We have 
proved that this curve is the set of singular points of this 2-surfacc characterized by the condition 

Q's2 = L's2 = R" = Q ^ Xs = Xe{z, a). 

This surface can be constructed using slices of z = constant. Each slice exhibits the limits 
of Q and L values, showing that beyond Za{a) no physical orbits (77 > 0) exist. It means that if 
E < Ea = m^fz^ neither spherical nor non null eccentricity orbits exist. The larger value of Za is 



4) Studying the \J potential we have seen that non-null eccentricity orbits are placed in this 
formal three-space between the two sheets of stable and unstable spherical orbits and limited by 
the surface representing the equatorial plane. Then, we obtained a global representation of the 
whole set of orbits in the outer region of the Kerr space-time. This representation allows us to see 
whether or not a particle can achieve one specific orbit from another one with a continuous change 
of its parameter values. 

5) Using a property of U potential, we have obtained xn^nit (see f|93p ). that is, a maximum 
value for the x coordinate of an orbit (spherical or not) depending on its energy-mass relation, E. 

6) We have obtained the necessary and sufficient condition that the variations of the Q, L and E 
parameters must hold so that a spherical orbit remains spherical. Using this condition, we recover 
the result obtained in |10j in the Schwarzschild limit and small perturbations of a stable spherical 
orbit, and in [12] in Kerr space-time. 

7) Variations of the parameters due to gravitational radiation in the adiabatical approximation 
can be represented in this formal space as a curve going form one geodesic to another. Therefore 
we can check whether or not this path is possible according to the fact that the upper sheet of S 
(that of stable spherical orbits) is impenetrable while the lower one is not, as has been proved in 
section 6. This representation, as a formalization of analytical results, can be an important tool 
to check whether a transition from one orbit to another is valid or not. 

In particular the necessary condition to go from one stable spherical geodesic to an unstable 
one (or vice versa) following a way that only include spherical orbits is crossing the singular line 
of S. 

8 Appendix I 

In order to prove that 77^1 < in P, we can do the next changes of variable and parameter: a) 
z = 2xs^^^/{2xs^^^ - 1), b) Xs =ts+ x+ and c) Xs^^^ = t^mux + x+ (then < t < t^maa;) where 
2^smax is defined in (j33p and x^ is the outer horizon. After applying this changes, we obtain 



3 for a = 1. 




^[J(a;,, a, z) + 2AV2a;,(2(l - z)xs + z)], 



(97) 



{2xs - lya- 



where 



J{xs,a, z) 



j{ts,a,t. 



■smax 
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f[2{tsmax - ts) + 4]tl + 



+ 4[6Vr^(t. -ts™a.) + (3(v/T^-2))]t2 



4[(1 - a^) + 4(1 - a^)ts,nax + Vl - a^]ts + 



a') 



(98) 



This implies that ry^i < in the intervals of interest, q.e.d. 



9 Appendix II 

After some manipulations we can see that the system of equations — Q's2 = is equivalent to 

= + G{x,,z), (99) 

where 

Ax 

G{xs,z) = ^{8{z-l)iz-2)xl + 3z{A^3z)xs + 3z^ + 



+ 4[2{z - l)xs - z]y/2xs[2{l - z)xs + z]}. (100) 
In fact, functions ^^.d 2^2 can be written 

e2(x.,a,z) = :^Ii^[a2 + G(x„z)], (101) 
a 

Q',2{xs,a,z) = ^^^^[a^ + G{xs,z)], (102) 

where F{xs, z) ^ and H{xs, z) 7^ 0. 

We define functions S{x, a, z) according to 

S{x,a,z) = G{x,z) + a^= (103) 



1^ - ^[8(^ - l){z - 2)x^ + 3z(4 - 2,z)x + 2,z\ 
16a; , 



- -^\2{z - l)x - z]J-2x\2{z - l)x - z\. (104) 
We implement the change of variable 

x^y = ^— -1/2, (105) 

where Xs,^^^ has been defined in ([33]). 

This change implies that .t = y = —1/2 and x = a;^^^^ ^ y = +1/2. At the same time 
y = ^ a; Xs^^^/2. In addition to this 

x+^y+ = l-- + il--)Vl-a^- (106) 
2 z z 

Now, the position of the outer horizon ?/-|_ depends not only on a but also on z. Remember that 
the domain V exists iff. condition x-^. 

<a;s„.ax: ^ 1 - (2 - 1)\/1 - a2 > holds. 
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The dependence of function S on y is 



S{y,z,a) = a^- \ [4{z - 2)y^ - 3zy^ + - 1] + (107) 
(z — 1)^ 4 

(l_4y2)3/2_ (^Qg) 



(Z- 1)3/2 

The values of S in the following points are 



5(y = -i) = a2>0, (109) 



^(. = +i) ^ a^-l + -^^^i^^^^^l^^>0, (110) 
^ 2^ (2-1)2 2(z- 1)^/1 ~ - ^ ' 



2 , z[5(2-l)-4(z-l)l/2 + l)] 



S{y^O) = -a' + ^ ' < ^^^^^ 

S{y+) = - - 1)^1 - + v/l - «2) X 

Z'^ 



X 1 - (z - l)\/l - a2 - Vl - a2]2} < 0. (112) 

As a direct consequence we see that S'(±l/2) > 0, S'(O) is a negative function of z, (1 < z), such 
that S'(O) < a2 - 1, and S'(y+) > 0. 
Moreover, 



6z?/ 

2-1 

The equation 5'(y) = has two solutions in 



^'(y) = r^[2(^ - 2)y - z + 2(z - l)V2(i _ 4y2)i/2]. (ug) 

(Z- 1)2 



yi = 0, 

V2 = 5(y2) = a2-l. (114) 



It is crucial to see that 



+ (l-l)(l-a2)V2 ^y+>y2, (115) 



and 5(0) = 5(2/i) < 5(^2) < 0. 

Considering these previous results, and taking into account the continuity of S", we can assert 
that this curve has at least one root [S = 0) in the interval — 1/2<?/<0 and at least another one 
in < y < +1/2. This is the proof of the existence of at least two solutions of the equation ([99]) . 
But we don't know if there exist solutions in the interval y+ < y < 1/2, and if there exist only one. 

The "mean- value theorem" states that if there exist two points ya,yb, Ha < 2/6, such that 
S{ya,b) = 0, then there is, at least, a point yo, ya < Uo < Ub, such that S'{yo) = 0. As we have at 
least two roots of the equation S{y, z, a) = 0, then, if we find yo, we can assert that one root must 
be a value greater than j/g (and another one less than yg). Therefore, if y+ < ?/o or y+ ^ Vi then 
we will prove that one single cut point will be outside the hole. 

We assert that there is a solution beyond the outer horizon iff. ?/-|_ < t/oi = 0, that is equivalent 

to 

1 + (1 - a2)i/2 

Z < ZUrmt(a) = (116) 

2 + (1 - a2)V2 
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since, according to (|115[) j/02 < y+- Nevertheless, there exist spherical orbits with z > zumit, so we 
don't know if in this limiting cases the root of S could be inside the outer horizon. But this is not 
possible since S{y+) < and then the root must be outside the horizon (the equality meets in the 
limit when j/+ — 1/2). 
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